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ABSTRACT 

We present a catalogue of accurate positions and correlated flux densities for 410 flat-spectrum, 
compact extragalactic radio sources previously detected in the Australia Telescope 20 GHz 
(AT20G) survey. The catalogue spans the declination range [—90°, — 40°] and was constructed 
from four 24-h very long baseline interferometry (VLB I) observing sessions with the Aus- 
tralian Long Baseline Array at 8.3 GHz. The VLBI detection rate in these experiments is 97 per 
cent, the median uncertainty of the source positions is 2.6 mas and the median correlated flux 
density on projected baselines longer than 1000 km is 0. 14 Jy. The goals of this work are (1) to 
provide a pool of southern sources with positions accurate to a few milliarcsec, which can be 
used for phase-referencing observations, geodetic VLBI and space navigation; (2) to extend 
the complete flux-limited sample of compact extragalactic sources to the Southern hemisphere; 
and (3) to investigate the parsec- scale properties of high-frequency selected sources from the 
AT20G survey. As a result of this VLBI campaign, the number of compact radio sources south 
of declination —40° which have measured VLBI correlated flux densities and positions known 
to milliarcsec accuracy has increased by a factor of 3.5. The catalogue and supporting material 
is available at http://astrogeo.org/lcsl. 
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1 INTRODUCTION 

Catalogues of positions of compact extragalactic radio sources with 
the highest accuracy are important for many applications. These 
include imaging faint radio sources in the phase-referencing mode, 
differential astrometry, space geodesy and space navigation. The 
method of very long baseline interferometry (VLBI) first proposed 
by Matveenko, Kardashev & Sholomitskii (1965) allows us to de- 
rive the position of sources with nanoradian precision (1 nrad ~ 
0.2 mas). The first catalogue of source coordinates determined with 
VLBI contained 35 objects (Cohen & Shaffer 1971). Since then 
hundreds of sources have been observed under geodesy and as- 
trometry VLBI observing programmes at 8.6 and 2.3 GHz (X and S 
bands) using the Mark3 recording system at the International VLBI 
Service for Geodesy and Astrometry (IVS) network. The analysis 
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of these observations resulted in the ICRF catalogue of 608 sources 
(Ma et al. 1998). 

The Very Long Baseline Array (VLB A) was later used to measure 
the positions of over 5000 compact radio sources in the VLBA 
Calibrator Survey (VCS) (Beasley et al. 2002; Fomalont et al. 2003; 
Petrov et al. 2005, 2006; Kovalev et al. 2007; Petrov et al. 2007) and 
the geodetic programme RDV (Petrov et al. 2009b). All sources with 
declinations above —45° detected using Mark3/Mark4 under IVS 
programmes were re-observed with the VLBA in the VCS and RDV 
programmes, which significantly improved the accuracy of their 
positions. As a result of these efforts, the probability of finding a 
calibrator with a VLBI-determined position greatly increased. In the 
declination range 8 > —30°, the probability of finding a calibrator 
within a 3° radius of a given position reached 97 per cent by 2008. 

Since the VLBA is located in the Northern hemisphere, obser- 
vations in the declination zone [—50°, — 30° ] are difficult and the 
array cannot observe sources with 8 < —52°. In 2008, the prob- 
ability of finding a calibrator within a radius of 3° was 75 per 
cent in the declination zone [—40°, — 30° ] and 42 per cent for 
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declinations south of —40°. The VLBI calibrator list 1 in 2008 had 
524 sources in the zone [+52°, +90°], but only 98 objects in the 
zone [—52°, — 90° ] which cannot be reached by the YLBA. These 
southern sources were observed during geodetic experiments and 
during two dedicated Southern hemisphere astrometry campaigns 
(Fey et al. 2004, 2006). The reason for this disparity is the scarcity 
of VLBI antennas in the Southern hemisphere, particularly stations 
with dual-frequency SIX receivers and geodetic recording systems. 
Also until recently there has been a lack of good all- sky catalogues 
suitable for finding candidate VLBI calibrators. 

The Australian Long Baseline Array (LBA) consists of six anten- 
nas located in Australia with the South Africa station HARTRAO 
often joining in. This VLBI network operates three to four observ- 
ing sessions a year, each about 1 week long. Although the hardware 
used by the LBA was not designed for geodesy and absolute astrom- 
etry observations, it was demonstrated by Petrov et al. (2009a) that 
despite significant technical challenges, absolute astrometry VLBI 
observations with the LBA network is feasible. In a pilot experiment 
in 2007 June, the positions of participating stations were determined 
with accuracies 3-30 mm and positions of five new sources were 
determined with accuracies 2-5 mas. 

Inspired by these astrometric results, we launched the X-band 
LBA Calibrator Survey (LCS) observing campaign, with the aim of 
determining milliarcsecond positions and correlated flux densities 
for 1000 compact extragalactic radio sources at declinations south of 
—30°. The overall objective of this campaign is to match the density 
of calibrators in the Northern hemisphere and so to eliminate the 
disparity. 

We have three long-term goals in this campaign. First, setting up a 
dense grid of calibrators with precisely known positions within sev- 
eral degrees of any target will make phase-referencing observations 
of weak targets feasible. According to Wrobel (2009), 63 per cent of 
VLBA observations in 2003-08 were made in the phase-referencing 
mode. A dense grid of calibrators also makes it possible to do dif- 
ferential astrometry of Galactic objects such as pulsars and masers, 
and allows the direct determination of the parallax at distances up to 
several kiloparsec (Deller et al. 2009). These sources form the pools 
of targets for observations under the geodesy programmes and for 
space navigation. 

Our second goal is to extend the complete flux-limited sam- 
ple of compact extragalactic radio sources (with emission from 
milliarcsecond- size regions) to the entire sky. According to Kovalev 
(private communication), an analysis of the log V-log S diagram of 
the VLBI calibrator list suggests that the sample of radio-loud ac- 
tive galactic nuclei (AGN) is complete at the level of the correlated 
flux density 200 mJy at X band at spatial frequencies 25 MX at 8 > 
—30°. Extending this complete sample to the entire sky will make 
it possible to generalize the properties of the sample, such as the 
distribution of compactness, distribution of the brightness temper- 
ature, bulk motion, viewing angle and irregularities of the spatial 
distribution, to the entire population of radio-loud AGN. 

The third goal is to investigate the properties of high-frequency 
selected radio sources from the Australia Telescope 20 GHz 
(AT20G) survey (Massardi et al. 2011; Murphy et al. 2010). Ob- 
taining observations of a subsample of AT20G sources with mil- 
liarcsecond resolution will allow us to investigate the properties, 
such as the spectral index, polarization fraction and variability, of a 
population of extremely compact sources. 


1 Available at http://astrogeo.org/rfc 


In this paper, we present the results from the first four 24-h ex- 
periments observed in 2008-09. The selection of candidate sources 
from the AT20G catalogue is discussed in Section 2. The station 
setups during the observing sessions is described in Section 3. The 
correlation and post-correlation analysis, which is rather different 
from ordinary VLBI experiments, is discussed in Sections 4 and 5. 
An error analysis of single-band observations, including the evalua- 
tion of ionosphere-driven systematic errors, is given in Section 5.2. 
The catalogue of source positions and correlated flux densities is 
presented in Section 6 and the results are summarized in Section 7. 

2 CANDIDATE SOURCE SELECTION 

The AT20G survey is a blind radio survey carried out at 20 GHz with 
the Australia Telescope Compact Array (ATCA) between 2004 and 
2008 (Murphy et al. 2010). It covers the whole sky south of decli- 
nation 0°. The source catalogue is an order of magnitude larger 
than previous catalogues of high-frequency radio sources, with 
5890 sources above a 20-GHz flux-density limit of 40mJy. All 
AT20G sources have the total intensity and polarization measured 
at 20 GHz, and most sources south of declination —15° also have 
near-simultaneous flux-density measurements at 4.8 and 8.6 GHz. 
A total of 1559 sources were detected in polarized total intensity at 
one or more of the three frequencies. There are also optical identi- 
fications and redshifts for a significant fraction of the catalogue. 

This high-frequency catalogue provides a good starting point 
for selecting bright, compact sources and candidate calibrators. 
Massardi et al. (201 1) show that almost all the flat-spectrum AT20G 
sources with spectral index a > —0.5 are unresolved on scales of 
0. 1-0.2 arcsec in size at 20 GHz (see their fig. 3.4). The few excep- 
tions are either (i) foreground Galactic and large magellanic cloud 
(LMC) sources such as planetary nebulae, H n regions and pulsar 
wind nebulae, which have a flat radio spectrum due to their thermal 
emission but are usually resolved on scales of a few arcsec (Mur- 
phy et al. 2010); or (ii) flat- spectrum extragalactic sources which 
are gravitationally lensed, like PKS 1830—211. 

For 64 per cent of sources from the AT20G catalogue, flux den- 
sities were determined in three bands, 5.0, 8.4 and 20 GHz. We 
used these measurements to calculate the spectral index a (F oc 
/“). We selected a set of 684 objects, not previously observed with 
VLBI, which had a flux density >150 mJy at 8.3 GHz and a spectral 
index > —0.6. We then split this sample into two subsets: 410 high- 
priority sources with flux densities >200mJy and spectral indices 
> —0.5, and all others. In addition, we selected 14 flat-spectrum 
objects that did not have flux-density measurements at 5 and 8 
GHz in the AT20G survey. Their spectral indices were determined 
by analysing historical single-dish observations found in the As- 
trophysical CATalogs support System (CATS, Verkhodanov et al. 
1997) data base, which by 2010 included data from 395 catalogues 
from radioastronomy surveys. 

In addition to these target sources, we identified a set of 195 
sources that we used for calibrator selection. These were bright 
sources previously observed at the VLBA or IVS network with 
positions known with accuracies better than 0.5 mas and with cor- 
related flux densities >500 mJy on baselines longer than 5000 km. 

2.1 The Australian Long Baseline Array 

The LBA network used for the LCS 1 consists of six stations: ATCA- 
104, CEDUNA, DSS45, HOBART26, MOPRA and PARKES 
(Fig. 1, Table 1). The maximum baseline length of the network 
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Figure 1 . The LBA network used for LCS1 observations. 

Table 1 . The LBA network stations. The average system equivalent flux 
density (SEFD) at 8.3 GHz at elevation angles >60° achieved in LCS1 
experiments is shown in the last column. 


Name 

0 gc 

k 

Diameter 

SEFD 

ATCA-104 

-30? 15 

149? 57 

22 m 

350 Jy 

CEDUNA 

— 31?70 

133?81 

32 m 

490 Jy 

DSS45 

— 35?22 

148?98 

34 m 

120 Jy 

HOBART26 

—42? 62 

147? 44 

26 m 

620 Jy 

MOPRA 

— 31? 10 

149? 10 

22 m 

380 Jy 

PARKES 

— 32282 

148?26 

64 m 

36 Jy 


is 1703 km, the maximum equatorial baseline projection is 1501 
km and the maximum polar baseline projection is 963 km. 


2.2 Observation scheduling 

Observation schedules were prepared with the software program 
sur_sked. The scheduling goal was to observe each target source 
with all antennas of the array in three scans of 120 s each in the first 
three experiments and in four scans in the last experiment. After 
tracking a target source for 120 s, each antenna immediately slewed 
to the next object. The minimum time-interval between consecutive 
scans of the same source was 3 h. The scheduling software computed 
the elevation and azimuth of each candidate source, with constraints 
on the horizon mask and maximum elevations to which antennas can 
point. It calculated the slewing time taking into account the antenna 
slewing rate, slewing acceleration and cable wrap constraints. The 
sources were scheduled to minimize the slewing time and fit the 
scheduling constraints. A score was assigned to each target source 
visible at a given time, depending on the slewing time, the history 
of past observations in the experiment and the amount of time that 
the source would remain visible to all antennas in the array. 

Every hour a set of four calibrator sources were observed: two 
scans where all antennas had an elevation in the range [12°, 45° ], 
one scan at elevations [32°, 45° ] and one scan at elevations [45°, 
90°]. Scans with calibrator sources were 70 s long. The scheduling 
algorithm for each set found all combinations of calibrator sources 
that fell in the right elevation range and selected the sequence of 


four objects that minimized the slewing time. The purpose of these 
observations was (1) to serve as amplitude and bandpass calibrators; 
(2) to improve the robustness of estimates of the path delay in the 
neutral atmosphere; and (3) to tie the positions of new sources to 
existing catalogues such as the ICRF catalogue (Ma et al. 1998). 

The scheduling algorithm assigned four scans for 18 per cent of 
target sources, three scans for 76 per cent of sources, two scans for 3 
per cent of sources and one scan for 3 per cent of sources. Of the 26 
targets with fewer than three scans in a given observing session, 18 
were observed in two experiments. During a 24-h experiment, 1 1.5— 
12.0 h was consumed in observing target sources and 1. 5-2.0 h in 
observing calibrators. 


3 OBSERVATIONS 

As mentioned in the introduction, not all of the telescopes in 
the LBA network are capable of observing in the typical geode- 
tic/astrometric mode of dual SIX frequency bands with multiple 
spaced sub-bands and recorded to a Mark5 VLB I system. 

The HOBART26, PARKES and DSS45 have a Mark5 recording 
system and a Mark4 (Whitney et al. 2004) baseband conversion 
rack. However, for the first epoch, the LBADR system (Phillips 
et al. 2009) was used. For all subsequent experiments, the Mark5 
system was used. Data were recorded on to normal Mark5 disk 
packs and shipped to the Max-Plank-Institut fur Radioastronomie 
in Bonn for processing. 

The ATCA-104, MOPRA and CEDUNA only have the standard 
LBA VLBI back-end consisting of an Australia Telescope National 
Facility (ATNF) Data Acquisition System (DAS) with an LBADR 
recorder. The ATNF DAS only allows two simultaneous intermedi- 
ate frequencies (IFs): either two frequencies or two polarizations. 
For each of these IFs, the input 64-MHz analog IF is digitally fil- 
tered to two contiguous 16-MHz bands. The ATCA-104, MOPRA 
and PARKES have two ATNF DASs; however, the IF conversion 
system at each telescope means it is impractical to run in any mode 
other than two frequencies and dual polarization. The LBADR data 
format is not compatible with the Mark4 data processor at Bonn. A 
custom program was written to translate the data to the Mark5B for- 
mat and then it was electronically copied using the Tsunami-UDP 
application to Bonn, before being copied on to Mark5 disk packs. 

In the last experiment, the PARKES recorded on the LBADR sys- 
tem in parallel with the Mark5 system. As this included left circular 
polarization (LCP) which was not recorded on the Mark5 system 
and since the ATCA-104 and MOPRA recorded both RCP and LCP, 
the LBADR data were also sent to Bonn and LCP correlated against 
the ATCA-104 and MOPRA. 

The main limiting factor for frequency selection was the LBADR 
back-end at the ATCA-104, MOPRA and CEDUNA where a setup 
with bands centred on ~8.2 and ~8.5 GHz was chosen. For the tele- 
scopes with Mark5 recorders, the setup was chosen to overlap in 
frequency with the LBADR setup but also including more frequen- 
cies to gain sensitivity. The frequency setup was changed between 
experiments in order to explore the feasibility of recording at 512 
Mbits -1 at those stations that can support it. The setup for each 
observing experiment is described in Table 2. 

The ATCA consists of six 22-m antennas and may observe as 
a single antenna or as a phased array of five dishes. 2 Since no 
tests of using a phased ATCA as an element of the VLBI network 


2 The sixth dish is located at the fixed pad far away from other antennas, 
which makes its phasing with the rest of the compact array too difficult. 
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Table 2. Observation summary for each observing session. 


Telescope Recorder Polarization Frequency bands 

(MHz) 


Experiment v254b, 2008 February 5 

PARKES 

ATCA-104 


MOPRA 

HOBART26 

CEDUNA 

LB ADR 

RCP only 

8256-8272, 8272-8288 
8512-8528, 8528-8544 

Experiment v271a, 2008 August 10 
PARKES 

HOBART26 Mark5 RCP only 

DSS45 

8200-8216, 8216-8232 
8232-8248, 8248-8264 
8456-8472, 8472-8488 
8488-8504, 8504-8520 

CEDUNA 

LB ADR 

RCP only 

8200-8216, 8216-8232 
8456-8472, 8472-8488 

ATCA-104 

MOPRA 

LB ADR 

RCP and LCP 

8200-8216, 8216-8232 
8456-8472, 8472-8488 

Experiment v271b, 2008 November 28 
PARKES 

HOBART26 Mark5 RCP only 

DSS45 

8200-8216, 8216-8232 
8264-8280, 8328-8344 
8392-8408, 8456-8472 
8472-8488, 8520-8536 

CEDUNA 

LB ADR 

RCP only 

8200-8216, 8216-8232 
8456-8472, 8472-8488 

ATCA-104 

MOPRA 

LB ADR 

RCP and LCP 

8200-8216, 8216-8232 
8456-8472, 8472-8488 

Experiment v271c, 2009 July 4 
HOBART26 Mark5 

RCP only 

8200-8216, 8216-8232 
8232-8248, 8328-8344 
8344-8360, 8456-8472 
8472-8488, 8488-8504 

PARKES 

Mark5 and 
LB ADR 

RCP and LCP 

8200-8216, 8216-8232 
8232-8248, 8328-8344 
8344-8360, 8456-8472 
8472-8488, 8488-8504 

CEDUNA 

LB ADR 

RCP only 

8200-8216, 8216-8232 
8456-8472, 8472-8488 

ATCA-104 

MOPRA 

LB ADR 

RCP and LCP 

8200-8216, 8216-8232 
8456-8472, 8472-8488 


for absolute astrometry observations were made before 2010, we 
decided to use a single dish of the ATCA in order to avoid the risk 
of introducing unknown systematic errors. 

The NASA Deep Space Network Station DSS45 observed only 
4.5 h in v271a and 6.5 h in v271b during intermissions between 
receiving the telemetry from the Mars-orbiting spacecraft. 

4 CORRELATION 

The Bonn Mark4 Correlator was chosen for the correlation of these 
experiments for two primary reasons. First, the correlator was ex- 
tensively tested for use in space geodesy and absolute astrometry 
mode during 2000-2010 and is known to produce highly reliable 
group delays. Secondly, it was equipped with four Mark5B units and 
eight Mark5 A units, which was a convenient combination since data 


from the PARKES were recorded in the Mark5B format and data 
from the ATCA-104, MOPRA and CEDUNA originally recorded in 
the LB ADR format were transformed to the Mark5B format before 
correlation. 

The preparation for the correlation required about 2 months due 
to complications that arose from the scheduling of different patch- 
ing and channel outputs (fan out) for use at the stations. The stations 
that had Mark4 data acquisition racks delivered detailed log files, 
from which we could reconstruct the channels used, but the LBA 
stations did not provide log files and the track assignments had to 
be searched for by trial and error, which was a time-consuming 
process. The absence of log files also required some custom pro- 
gramming to reconstruct the log VEX file (lvex), which the Mark IV 
correlator required to perform the correlation. Further, the setup that 
was chosen was incompatible with the capabilities of the hardware 
correlator and required an extensive workaround at the correlator. 

We chose to correlate with a window of 128 lags (corresponding 
to a delay window width of 8 |xs) instead of the 32 lags (2 |xs) nor- 
mally used for stream correlation. This was to allow for potentially 
large clock offsets due to instrumental errors and for large a priori 
source coordinate errors. The integration time was chosen to be as 
small as possible (0.5-1. 0 s) to allow for potentially large residual 
fringe rates due to source position errors. 

The Mark4 correlator can cope with only four stations simultane- 
ously when correlating 128 lags with short integration time. Since 
the experiment was observed with five to seven stations in right 
circular polarization (RCP) and also two stations (ATCA-104 and 
MOPRA) in LCP, the correlation had to be split in passes. 

As a by-product of the multipass correlation, some baselines 
were correlated more than once and this redundancy was exploited 
to select the correlation with the best signal-to-noise ratio (S/N) for 
each duplicated scan. For the LCP correlation, only one pass was 
required since only the LBA stations recorded LCP. The PARKES 
recorded LCP only during v271c and this was enabled by the use 
of two back-ends (Mark4 and LBA) in parallel. Those stations 
were recorrelated in a second pass to produce all four polarization 
products (left-left, right-right, left-right and right-left). For this 
pass, we had to restrict the number of lags to 32 due to correlator 
hardware constraints, but the integration time was kept at 0.5 s. 

Fringe fitting was performed at the correlator using the soft- 
ware program fourfit, the baseline-based fringe fit offered within 
the hops (Haystack Observatory Package Software) to estimate the 
residual delay upon which the post-correlation analysis was based. 
Inspection of fringe-fitted data is a convenient means of data-quality 
control, as one can examine the correlated data on a scan-by-scan 
basis, looking for problems that might have occurred during correla- 
tion or at the stations and permits flagging of bad data immediately. 


5 DATA ANALYSIS 

The correlator either generates the spectrum of the cross-correlation 
function directly or can be easily derived from its raw output. 
The spectrum was processed with the software program fourfit, 
which for each scan and each baseline determined the phase-delay 
rate, narrow-band delay and wide-band group delay (sometimes 
also called multiband delay) that maximized the fringe amplitude. 
Takahashi et al. (2000) give a detailed description of the fringe 
searching process and the distinction between narrow-band and 
wide-band group delays. The wide-band delay is more precise than 
narrow-band delay. The formal uncertainties of these delays, o n and 
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Table 3. The list of target sources that were not detected in LCS 1 VLBI observations. The last three data columns contain 
the flux density in mJy extrapolated to 8.3 GHz (F), the spectral index (a) and the comment (C). Value —9.99 means 
that the spectral index was not available. Comments: (1) not in the AT20G catalogue; (2) flagged as an LMC source in 
the AT20G catalogue; (3) Galactic planetary nebula; and (4) flagged as an extended source in the AT20G catalogue. This 
table is also available electronically (see Supporting Information). 


J2000 name 

B1950 

RA 

Dec. 

F 

a 

C 

J0404— 7109 

0404-712 

04 04 00.99 

-71 09 09.7 

111 

-0.57 

(1) 

J0538— 6905 

0539-691 

05 38 45.66 

-69 05 03.1 

204 

-9.99 

(2) 

J0552— 5349 

0551-538 

05 52 36.18 

-53 49 32.4 

182 

-0.34 

(1) 

J0938— 6005 

0937-598 

09 38 47.20 

-60 05 28.7 

196 

-0.08 

(3) 

J0958— 5757 

0956-577 

09 58 02.92 

-57 57 42.6 

383 

0.21 

(3) 

J1100— 6514 

1058-649 

11 00 20.09 

-65 14 56.4 

179 

-0.08 

(3) 

J1150— 5710 

1147-569 

11 50 17.87 

-57 10 56.0 

349 

0.23 

(3) 

J1325— 4302 

1322-427 

13 25 07.35 

-43 02 01.8 

327 

-9.99 

(4) 

J1353— 6630 

1350-662 

13 53 57.05 

-66 30 50.3 

370 

0.04 

(3) 

J1505— 5559 

1502-557 

15 05 59.17 

-55 59 16.2 

111 

-0.06 

(3) 

J1656— 4014 

1653-401 

16 56 47.53 

-40 14 24.4 

All 

-9.99 

(4) 


o W9 are computed by fourfit in the following way: 

Vl2 

° n ~ 2n A/ (S/N) ’ 

a* = , (1) 

2nof (S/N) 

where A/ is the IF bandwidth, Of is the dispersion of cyclic fre- 
quencies across the band and S/N is the signal-to-noise ratio of the 
fringe amplitude from the wide-band fringe search to the rms of 
the thermal noise. The ratio of a Jo n was in the range 27-28 for 
the LCS1 experiments, which gives for observations with typical 
S/N = 30 <j w ^ 40 ps and o n around 1 ns. 

Of the 421 observed targets, 410 were detected. The list of 11 
undetected sources is given in Table 3. In addition, 111 calibrators 
were observed, all of which were detected. 

5.1 Data analysis: source position determination 

The most challenging part of the data analysis was resolving group- 
delay ambiguities. The algorithm for fringe fitting implemented in 
fourfit searches for a global maximum in the Fourier transform 
of the cross-correlated function, averaged over individual IFs after 
correcting phases for a fringe-delay rate and a narrow-band group 
delay. The fringe spectrum folding results in a rail of maxima of 
the Fourier transform of the cross-correlated function with exactly 
the same amplitude and with spacings reciprocal to the minimum 
difference in intermediate frequencies, 1/16.0 x 10 6 Hz = 62.5 ns. 
Within one half of the 62.5 ns range, the fringe spectrum has sev- 
eral strong secondary maxima. The second maximum at 1/256.0 x 
10 8 Hz & 3.9 ns has the amplitude 0.981 of the main maximum and 
the third maximum at 2/256.0 x 10 8 Hz ^ 7.8 ns has the amplitude 
0.924. This high level of secondary maximum amplitudes is due to 
our choice of IFs that was determined by the hardware limitations. 
Due to the presence of noise in the cross-correlated function, the 
difference between the fringe amplitude at the main maximum and 
at the secondary maxima is random and therefore the group delay 
is determined with an uncertainty ±1-3 of the spacing between the 
secondary maxima, 3.9 ns. 

It should be noted that typical group-delay ambiguity spacings in 
geodetic observations are in the range 50-200 ns. For successfully 
resolving ambiguities, as a rule of thumb, the predicted delay should 
be known with an accuracy better than one- sixth of the ambiguity 
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Figure 2. The postfit residuals of the narrow-band group-delay LSQ solu- 
tion at the baseline CEDUNA/HOBART26 in experiment v271a on 2008 
February 5. The bar in the middle of the plot corresponds to the ambiguity 
spacing of the wide-band path group delay. 

spacing, 600 ps in our case. This is a challenge. The random errors 
of the narrow-band group delays are too high to use directly for 
resolving ambiguities in wide-band group delays (see Fig. 2). In 
the framework of the traditional data analysis of geodetic VLBI 
observations, the maximum uncertainty in the prediction of the 
path delay is due to a lack of an adequate model for the path delay 
in the neutral atmosphere. These errors are in the range of 300- 
2000 ps. The path delay through the ionosphere at 8.3 GHz is in the 
range 30-1000 ps. In addition, an error ofl arcsec in a priori source 
position (a typical AT20G position error) causes an error in a priori 
time-delay of up to 30 ns on a 1700 km baseline. 


5.1.1 Group-delay ambiguity resolution 

However, it is premature to conclude that resolving group-delay 
ambiguities is impossible. First, we need to use a state-of-the-art 
a priori model. Our computation of theoretical time-delays in gen- 
eral follows the approach of Sovers, Fanselow & Jacobs (1998) 
with some refinements. The most significant ones are the following. 
The advanced expression for the time-delay derived by Kopeikin & 
Schafer (1999) in the framework of general relativity was used. The 
displacements caused by the Earth’s tides were computed using 
the numerical values of the generalized Love numbers presented 
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by Mathews (2001) following a rigorous algorithm described by 
Petrov & Ma (2003) with a truncation at a level of 0.05 mm. The 
displacements caused by ocean loading were computed by convolv- 
ing Greens’ functions with ocean tide models. The GOT99.2 model 
of diurnal and semi-diurnal ocean tides (Ray 1999), the NA099 
model (Matsumoto, Takanezawa & Ooe 2000) of ocean zonal tides, 
the equilibrium model (Petrov & Ma 2003) of the pole tide and 
the tide with a period of 18.6 yr were used. Station displacements 
caused by the atmospheric pressure loading were computed by con- 
volving Greens’ functions that describe the elastic properties of the 
Earth (Farrell 1972) with the output of the atmosphere NCEP Re- 
analysis numerical model (Kalnay et al. 1996). The algorithm of 
computations is described in full detail in Petrov & Boy (2004). 
The displacements due to loading caused by variations in soil mois- 
ture and snow cover in accordance with the GLDAS Noah model 
(Rodell et al. 2004) with a resolution 0.25 x 0.25 deg 2 were com- 
puted using the same technique as the atmospheric pressure loading. 
The empirical model of harmonic variations in the Earth orienta- 
tion parameters heo_20101111 derived from VLBI observations 
according to the method proposed by Petrov (2007) was used. The 
time-series of UT1 and polar motion derived by the NASA God- 
dard Space Flight Center operational VLBI solutions were used a 
priori. 

The a priori path delays in the neutral atmosphere in the direc- 
tion of observed sources were computed by the numerical inte- 
gration of differential equations of wave propagation through the 
heterogeneous media. The four-dimensional field of the refractivity 
index distribution was computed using the atmospheric pressure, 
air temperature and specific humidity taken from the output of the 
Modern Era Retrospective- Analysis for Research and Applications 
(MERRA) (Schubert et al. 2008). That model presents the atmo- 
spheric parameters at a grid 1/2° x 2/3° x 6 h at 72 pressure levels. 

Secondly, we made a least-squares (LSQ) solution using narrow- 
band delay. The positions of all target sources were estimated, as 
well as clock functions for all stations except the one taken as a ref- 
erence and the residual atmosphere path delay in zenith directions. 
The clock function was modelled as a sum of the second-degree 
polynomial over the experiment and B- spline of the first order with 
the time-span 60 min. The residual atmosphere path delay in the 
zenith direction was modelled with B -spline of the first order with 
a time-span 60 min. Constraints on the rate of change of the clock 
function and atmosphere path delay were imposed. After the re- 
moval of outliers (1-2 per cent of points), the weighted rms (wrms) 
of residuals was 1-4 ns. An example of narrow-band postfit resid- 
uals is shown in Fig. 2. 

The next step was to update the differences between the wide- 
band path delays and the theoretical group delays by subtracting the 
contribution to path delay of the adjustments from the narrow band 
LSQ solution to clock function, atmosphere path delay in zenith 
direction and source positions. An example of enhanced wide-band 
group delays after the substitution and group delay ambiguity res- 
olution is shown in Fig. 3. 

It is still difficult to resolve ambiguities at long baselines, but 
relatively easy to resolve the ambiguities at the inner part of the ar- 
ray: ATCA-104, DSS45, MOPRA and PARKES. The group-delay 
ambiguity resolution process starts from the baseline with the least 
scatter of a priori wide-band delays. After the group-delay am- 
biguity resolution at the first baseline and temporary suppression 
of observations with questionable ambiguities, the LSQ solution 
with mixed delays is made: wide-band group delays at baselines 
with resolved ambiguities and narrow-band group delays at other 
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Figure 3. The wide-band a priori group delays after the substitution of 
adjustments to the clock function, atmosphere path delay in the zenith direc- 
tion and source positions from the narrow-band LSQ solution at the baseline 
CEDUNA/PARKES in experiment v271a on 2008 February 5. The bar on 
the left-hand side of the plot corresponds to the ambiguity spacing of the 
wide-band path group delay. 

baselines. In addition to other parameters, baseline-dependent clock 
misclosures are estimated. 

In the absence of instrumental delays, the differences between 
wide-band and narrow-band delays would be the zero mean Gaus- 
sian random noise. Instrumental delays in the analogue electronics 
cause systematic changes of these differences in time. Fortunately, 
these systematic changes are smooth and small enough to allow 
solving for ambiguities with spacings as small as 3.9 ns. Since the 
scatter of the postfit residual of wide-band group delays is a factor 
of 10-50 less than the ambiguity spacings, when the instrumental 
delay is determined, the ambiguities can be resolved. Therefore, 
the feasibility of resolving wide-band group delays hinges upon the 
accuracy of the determination of the differences between narrow- 
band and wide-band group delays from the analysis of narrow-band 
delays. 

An example of the a priori wide-band group delay with the sub- 
stituted adjustments to source positions found from the wide-band 
group delay solution at other baselines is shown in Fig. 4. After 
resolving ambiguities for the most of observations, the points that 
were temporarily suppressed were restored and used in the solution. 



Figure 4. The wide-band a priori group delays at the baseline CE- 
DUNA/PARKES in experiment v271a on 2008 February 5 after the sub- 
stitution of adjustments to the clock function, atmosphere path delay in the 
zenith direction from the narrow-band LSQ solution and source positions 
from the wide-band LSQ solution at other baselines and group-delay ambi- 
guity resolution. The bar in the middle part of the plot corresponds to the 
ambiguity spacing of the wide-band path group delay. 
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5.1.2 Ionosphere path-delay contribution 

Single-band VLB I data are affected by a variable path delay through 
the ionosphere. We attempted to model this path delay using GPS 
TEC maps provided by the CODE analysis center for processing 
Global Navigation Satellite System data (Schaer 1998). Details of 
computing the VLBI ionospheric path delay using TEC maps from 
the GPS are given in Petrov et al. (201 1). The model usually recovers 
over 80 per cent of the path delay at baselines longer than several 
thousand kilometres. However, applying the ionosphere TEC model 
to processing the LCS1 experiment did not improve the solution 
and even degraded the fit. The same model and software program 
certainly improved the fit and improved results when applied to 
processing observations on intercontinental baselines. 

In order to investigate the applicability of the reduction for the 
path delay in the ionosphere based on GPS TEC maps, we processed 
29 IVS dual-band geodetic experiments that included the 1089-km- 
long baseline HOBART26/PARKES. We ran three solutions that 
included the data only at this baseline. The first reference solution 
used ionosphere-free linear combinations of X-band and 5-band 
observables. The second solution used X-band-only group delays. 
The third solution used X-band-only data and applied the reduction 
for the path delay in the ionosphere from GPS TEC maps. 

We discarded two experiments: one had a clock break at HO- 
BART and another had the rms of a postfit residual a factor of 5 
greater than usually due to a warm receiver. The rms differences of 
postfit residuals with respect to the reference dual-band X/S solution 
for two trial solutions are shown in Fig. 5: one with X-band-only 
data and another with X-band data with ionosphere path delays from 
GPS TEC maps applied. Considering the reference solution based 
on ionosphere-free group delays as the ground truth, we expected 
the rms differences between the X-band-only solution (denoted with 
the circles in Fig. 5) to be positive and the rms differences of the 
X-band solution with the ionosphere path delay from GPS TEC 
maps applied to be also positive but less than the rms differences 
between the X-band-only solutions. Instead, we see that the rms of 
postfit residuals of the X-band-only solution is less than the rms of 
the dual-band reference solution and applying the ionosphere path 
delay from GPS TEC maps does not affect the rms significantly. 

We then computed the rms of the contribution of the ionosphere 
path delay from dual-band X/S observations and the rms of the 
differences between the path delay in the ionosphere from the dual- 
band X/S observations and the GPS TEC ionosphere maps. They are 



2002 2004 2006 2008 2010 

Time in years 

Figure 5. The difference in the rms of postfit residuals with respect to the 
reference dual-band X/S solution for two trial solutions: X-band-only data 
(solid discs) and X-band data with the ionosphere path delay from GPS TEC 
maps applied (circles). 



Figure 6. The rms of the ionosphere contribution from dual-band X/S ob- 
servations (solid discs) and the rms of the differences between the path 
delay in the ionosphere from the dual-band X/S observations at the baseline 
HOBART26/PARKES and the GPS TEC ionosphere maps (circles). 

shown in Fig. 6. We see that in 2002-2004 the ionosphere path delay 
from GPS TEC maps was coherent with the dual-band ionosphere 
path delay estimate. In 2005-2010, the ionosphere path delay was 
small with rms ~30 ps and the ionosphere path delay from the GPS 
TEC model was only partially coherent with the ionosphere path 
delay estimates from dual-band VLBI observations. The baseline 
length repeatability, defined as the rms of baseline length estimates 
after subtracting the secular drift due to tectonics, is the minimum 
when ionosphere-free linear combinations of X/S observables are 
used: 17.0 mm, grows to 17.4 mm when the X-band group delays 
are used and it is the maximum, 18.4 mm, when the ionosphere 
path-delay reduction based on the GPS TEC model is applied to the 
X-band observables. 

We compared these results with analogous observations in the 
Northern hemisphere. We picked the baseline LA-VLBA/OV- 
VLBA that has almost exactly the same length as the baseline 
HOBART26/PARKES: 1088 km against 1089 km. We processed 
88 experiments at this baseline in a similar fashion to the way we 
processed the HOBART26/PARKES baseline. Results are shown 
in Fig. 7. Applying the ionosphere path delay from the GPS has 
reduced the rms of postfit residuals at this baseline. 

Our interpretation of this phenomenon is that the ionosphere path 
delay from GPS TEC maps at baselines 1000-2000 km has the floor 
around 30 ps. The dominant constituents in the ionospheric path 



Figure 7. The rms of the ionosphere contribution from dual-band X/S ob- 
servations (solid discs) and the rms of the differences between the path delay 
in the ionosphere from the dual-band X/S observations at the 1088-km-long 
baseline LA-VLBA/OV-VLBA and the GPS TEC ionosphere maps (circles). 
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delay at these scales during the solar minimum are short-periodic 
scintillations that the GPS TEC model with time-resolution 2 h and 
spatial resolution 500 km does not adequately represent. It is also 
known (Hemandez-Pajares et al. 2009) that the accuracy of the GPS 
TEC model is worse in the Southern hemisphere than in the Northern 
hemisphere due to the disparity in the GPS station distribution. 

Applying the GPS TEC map ionosphere path-delay reduction 
effectively added noise in the data. Therefore, we did not apply the 
ionosphere path delay in our final solution. 


5.1.3 Reweighting observations 


According to the Gauss-Markov theorem, the estimate of param- 
eters has the minimum dispersion when observation weights are 
chosen reciprocal to the variance of errors. The group delays used 
in the analysis have errors due to the thermal noise in fringe phases 
and due to the mismodelling of the propagation delay: 


a 2 = a t l + <r n 2 a 


+ a in 


(2) 


where a th is the variance of the thermal noise, and cr na and o l0 are 
the variances of errors of modelling the path delay in the neutral 
atmosphere and the ionosphere, respectively. 

A rigorous analysis of the errors of modelling the path delay in 
the neutral atmosphere is beyond the scope of this paper. Assuming 
that the dominant error source of the a priori model is the high- 
frequency fluctuations of the water vapour at time- scales less than 
3-5 h, we sought a regression model for the dependence of a * a on 
the non-hydrostatic component of the slanted path delay through the 
atmosphere. We made several trial runs using all observing sessions 
under geodesy and absolute astrometry programmes for 30 yr and 
added in quadrature to the a priori uncertainties of group delay an 
additive correction: 


— + 



2 


( 3 ) 


where r s is the contribution of the non-hydrostatic constituent of the 
slanted path delay, r z is the non-hydro static path delay in the zenith 
direction computed by a direct integration of equations of wave 
propagation through the atmosphere using the refractivity computed 
using the MERRA model and a is a coefficient. We found that the 
baseline length repeatability defined as the rms of the deviation of 
the baseline length with respect to the linear time-evolution reaches 
the minimum when a = 0.02. Therefore, we adopted this value in 
our analysis of LCS1 experiments. For typical values of r z , the 
added noise is 8 ps in the zenith direction and 16 ps at the elevation 
of 30°. 

We computed the rms of ionosphere contributions from dual-band 
VLBI group delays and from GPS TEC maps. Fig. 8 shows the de- 
pendence of the square of the rms from dual-band VLB I versus the 
square of the rms from the GPS. The slope of the regression straight 
line is 0.992 =b 0.004. This dependence suggests that although the 
ionosphere path delay from TEC maps from the GPS analysis is 
too noisy to be applied for the reduction of observations at short 
baselines during the solar minimum, it correctly predicts the vari- 
ance of the ionosphere path delay. We assigned the variance of the 
mismodelled ionosphere path delay to the ionosphere contribution 
to the delay computed from TEC maps: o l0 = r io . 

We also computed, for each experiment and for each baseline, ad 
hoc variances of observables that after being added in quadrature 
make the ratio of the weighted sum of squares of postfit residuals 
to their mathematical expectation close to unity. This computation 
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Figure 8. The dependence of the ionosphere path delay cr? from dual-band 
VLBI versus from TEC maps from the GPS analysis. Data in the time 
range 2005.0-2010.6 (solar minimum) were used. 


technique is presented in Petrov et al. (2011). The ad hoc vari- 
ance was applied to further inflate the formal uncertainties of the 
observables that have already been corrected for the inaccuracy 
of the a priori model of wave propagation through the ionosphere 
and neutral atmosphere (equation 2). In contrast to cr aa and cr3, the 
baseline-dependent ad hoc variance is elevation-independent. 


5.1.4 Global LSQ solution 

We ran a single global LSQ solution using all available dual-band 
VLBI observations under geodesy and absolute astrometry pro- 
grams from 1980.04.12 through 2010.08.04, in total 7.56 million 
observations, and the X-band VLBI data from four LCS1 observing 
sessions. The RCP and LCP data were treated as independent exper- 
iments. The following parameters were estimated over the global 
data set: coordinates of 4924 sources, including 410 detected target 
objects in the LCS1 campaign (see Fig. 9); positions and velocities 
of all stations; coefficients of the expansion over the B- spline basis 
of non-linear motions for 17 stations; coefficients of harmonic site 
position variations of 48 stations at four frequencies: annual, semi- 
annual, diurnal and semidiurnal; and axis offsets for 67 stations. In 
addition, the following parameters were estimated for each exper- 
iment independently: station-dependent clock functions modelled 
by second-order polynomials, baseline-dependent clock offsets, the 
pole coordinates, UT1-TAI and daily nutation offset angles. The 
list of estimated parameters also contained more than one million 
nuisance parameters: coefficients of linear splines that model the at- 
mospheric path delay (20-min segment) and clock function (60-min 
segment). 

The rate of change for the atmospheric path delay and clock 
function between adjacent segments was constrained to zero with 
weights reciprocal to 1.1 x 10 -14 and 2 x 10 -14 , respectively, in 
order to stabilize our solution. We apply no-net rotation constraints 
on the position of 212 sources marked as ‘defining’ in the ICRF 
catalogue (Ma et al. 1998) that requires the positions of these sources 
in the new catalogue to have no rotation with respect to the position 
in the ICRF catalogue. 

The global solution sets the orientation of the array with respect 
to an ensemble of ~5000 extragalactic remote radio sources. The 
orientation of that ensemble is defined by the series of the Earth 
orientation parameters evaluated together with source coordinates. 
Common 111 sources observed in the LCS1 as atmosphere and 
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Figure 9. The distribution of sources with milliarc second positions at the 
southern sky. The circles denote sources known by 2008. The filled green 
discs denote sources from the LCS1 catalogue. The Galactic plane is shown 
by the red line. 

amplitude calibrators provide a strong connection between the new 
catalogue and the old catalogue of compact sources. 


5.2 Error analysis of the LCS1 catalogue 

To assess the systematic errors in our results, we exploited the 
fact that 111 known sources were observed as amplitude and at- 
mospheric calibrators. Positions of these sources were determined 
from previous dual-band SIX observations with accuracies better 
than 0.2 mas. We sorted the set of 1 1 1 calibrators according to their 
right ascensions and split them into two subsets of 55 and 56 objects, 
even and odd. We ran two additional solutions. In the first solution, 
we suppressed 55 calibrators in all experiments but the LCS1 and 
determined their positions solely from the LCS1 experiment. In the 
second solution, we did the same with the second subset. Consider- 
ing that the positions of calibrators from numerous SIX observations 
represent the ground truth, we treated the differences as LCS1 er- 
rors. 

We computed the x 2 /ndf statistics for the differences in right 
ascensions and declinations, A a and A s , and sought additional 
variances, v a and that, being added in quadrature to the for- 
mal source position uncertainties, a a j and cr s ,i, make them close to 
unity: 

xl = ^2 Act? cos 2 & 

ndf n 22 \ J a «.‘ cos2 

xf_ = 22 AlS - 2 (4) 

ndf n 22 + v * 

We found the following additive corrections of the uncertainties 
in right ascensions scaled by cos 8 and for declinations, v a = 1.44 
mas and v s = 0.51 mas, respectively. We do not have an explanation 
why the additive correction for scaled right ascensions is three times 


greater than for declinations. After applying the additive corrections, 
the wrms of source position differences is 1.8 mas for the right 
ascension scaled by cos 8 and 1.5 mas for declinations. 

The final inflated uncertainties of source positions, a 2 (/) and 
cr 2 (/), are, respectively, 

cr^f) = cr^ + v 2 a /cos 2 S 

and 

+ vl (5) 


5.3 Data analysis: correlated flux-density determination 

Each detected source has from 3 to 60 observations, with a median 
value of 25. Imaging a source with 25 points at the uv plane is a 
difficult problem and the dynamic range of such images will be 
between 1:10 and 1:100, which is far from spectacular. Images 
produced with the hybrid self-calibration method will be presented 
in a separate paper. In this study, we limited our analysis to mean 
correlated flux-density estimates in three ranges of lengths of the 
baseline projections on to the plane tangential to the source, without 
the inversion of calibrated visibility data. Information about the 
correlated flux density is needed for the evaluation of the required 
integration time when an object is used as a phase calibrator. 

First, we calibrated raw visibilities v for the a priori system 
temperature T sys and antenna gain G(e): A corr = vT sys (t , e)/G(e). The 
coefficients of antenna gain expansions into polynomials over the 
elevation angle e are presented in Table 4. The system temperature 
was measured at each station, each scan. 

In the second step, we adjusted antenna gains using publicly 
available brightness distributions of calibrator sources made with 
observations under other programmes. We compiled the data base 
of brightness distributions, correlated flux densities and images of 
compact radio sources produced with VLB I 3 from authors who 
agreed to make their imaging results publicly available. Among 
111 sources observed as calibrators, images of 14-27 objects at 
each experiment were available. These are images in the form of 
CLEAN components mainly from Ojha et al. (2004, 2005). CLEAN 
components of source brightness distributions from the analysis 
of the TANAMI program (Ojha et al. 2010), that observed with 
the LBA concurrently with the LCS1 were not available, but the 
parameters of one-component Gaussian models that fit the core 
regions were published. For those sources for which both a set 
of CLEAN components and the parameter of the Gaussian one- 
component model were available, we used CLEAN components. 

We predicted the correlated flux density F corr for each observation 
of an amplitude calibrator with known brightness distribution as 


Fron — 




- (ux+vy) 


p 4 TFT [ fl2 ( u cosfl+D sin0) 2 +b 2 (— u sin0+u cos#) 2 ] 


( 6 ) 


where c t is the correlated flux density of the ith CLEAN component 
with coordinates x and y with respect to the centre of the image, u 
and v are the projections of the baseline vectors on the tangential 
plane of the source, a and b are the full width at half-maximum of 
the Gaussian that approximates the core, and 0 is the position angle 
of the semimajor axis of the Gaussian model. 


3 Available at http://astrogeo.org/vlbi_images 
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Table 4. The coefficients of antenna gain polynomials over the elevation angle expressed in radians. Gain = DFPU J2k= o a k ek - The degrees per unit flux 
(DFPU) is shown in the second column. Columns (3)-(8) hold coefficients of degree 0-5. 


Station DFPU (K Jy l ) The coefficients of the polynomial for gain as a function of the elevation angle in degrees 

0 1 2 3 4 5 


ATCA-104 

0.100 

1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

CEDUNA 

1.000 

0.9045 

4.032 x 10“ 3 

-4.280 x 1(T 5 

0.0 

0.0 

0.0 

DSS45 

0.240 

0.9700 

2.927 x 1(T 3 

-1.070 x 1CT 4 

1.761 x 10 -6 

-1.112 x 10" 8 

0.0 

HOBART26 

0.058 

0.6997 

4.445 x 10" 2 

-2.331 x 1CT 3 

5.335 x 1CT 5 

-5.412 x 10- 7 

1.926 x 10“ 9 

MOPRA 

0.095 

1.0 

0.0 

0.0 

0.0 

0.0 

0.0 

PARKES 

0.480 

0.7900 

3.368 x 10“ 3 

8.507 x 10“ 5 

-1.398 x 1(T 6 

0.0 

0.0 


Table 5. The first 12 rows of the catalogue of source-correlated flux-density estimates for both 410 target sources and 111 calibrator sources in each experiment. 
Some sources were observed on more than one experiment. The table columns are explained in the text. The full table is available in the electronic version of 
this paper (see Supporting Information). 


J2000 name 

B1950 name 

Date 

Status 

# Observation 

Correlated flux density 

(Jy) 

Flux-density uncertainties 

(Jy) 

Exp 

J0004— 4736 

0002-478 

2008.11.28 


12 

0.352 

0.268 

0.267 

0.015 

0.019 

0.013 

v271b 

J001 1—8443 

0009-850 

2008.02.05 


30 

0.161 

0.165 

0.180 

0.005 

0.005 

0.006 

v254b 

J0012— 3954 

0010-401 

2008.02.05 


16 

0.761 

0.769 

0.715 

0.016 

0.012 

0.014 

v254b 

J0012— 3954 

0010-401 

2009.07.04 

Cal 

10 

0.815 

0.654 

0.640 

0.026 

0.023 

0.028 

v271c 

J0028— 7045 

0026-710 

2008.11.28 


40 

0.092 

0.086 

0.081 

0.005 

0.006 

0.005 

v271b 

J0030— 4224 

0027-426 

2009.07.04 


30 

0.312 

0.290 

0.256 

0.011 

0.010 

0.011 

v271c 

J0033— 4236 

0031-428 

2009.07.04 


30 

0.132 

0.135 

0.143 

0.005 

0.005 

0.008 

v271c 

J0034— 4116 

0031-415 

2009.07.04 


18 

0.053 

0.045 

0.039 

0.003 

0.002 

0.004 

v271c 

J0040— 4253 

0038-431 

2009.07.04 


28 

0.068 

0.058 

0.067 

0.003 

0.002 

0.007 

v271c 

J0042— 4030 

0039-407 

2009.07.04 


26 

0.255 

0.253 

0.276 

0.010 

0.011 

0.013 

v271c 

J0042— 4333 

0040-438 

2009.07.04 


30 

0.165 

0.113 

0.126 

0.005 

0.005 

0.008 

v271c 

J0044— 8422 

0044-846 

2008.02.05 


26 

0.278 

0.255 

0.234 

0.007 

0.006 

0.007 

v254b 


Then we built a system of LSQ equations for all observations of 
calibrators with known images used in astrometric solutions: 

^corr — V^I7 A corr (7) 

and after taking logarithms from the left- and right-hand sides, we 
solved for corrections to gains g for all stations. Finally, we applied 
corrections to gain for observations of all other sources. 

The correlated flux density is a constant during the observing 
session only for unresolved sources. For resolved sources, the cor- 
related flux density depends on the projection of the baseline vector 
on the source plane and on its orientation. We binned the corre- 
lated flux densities in three ranges of the baseline vector projection 
lengths, 0-6 MX, 6-25 MX and 25-50 MX, and found the median 
value within each bin. These bins correspond to scales of the de- 
tected emission at >30 mas, 7-30 mas and <7 mas, respectively. 

The uncertainties of our estimates of the correlated density de- 
pend on the thermal noise described as cr th = F corr /(S/N) and errors 
of gain calibration. The variance of the thermal noise was in the 
range of 1-6 mJy depending on the sensitivity of a baseline, with 
the median value of 3 mJy. The LSQ solution for gains provided the 
variance of the logarithm of gains. Assuming the calibration errors 
to be multiplicative in the form g (1 + e), where g is gain and e is 
the Gaussian random variable, we can evaluate the contribution of 
gain errors in the multiplicative uncertainty of F corr (l + /x), where 
/x the Gaussian random variable. Its variance is evaluated as 

cr(/x) = exp yj Cov(vi, si) + Covfe, s 2 ) + 2Cov(vi, s 2 )j , 

V J (8) 


where Co \{s t , sj) is the covariance of the logarithm of gain between 
the ith and yth stations. Multiplicative gain uncertainties are in the 
range 0.08-0.1 for the DSS45 and in the range 0.02-0.05 for other 
stations. The gain uncertainties for the DSS45 are higher, because 
it observed three to five times less than other stations. The total 
variance is a sum in quadrature of <r c and cr th . We should note 
that our estimate of systematic errors does not account for possible 
errors in the gain curve determination. A systematic error in the 
gain curve would directly affect our estimate of the correlated flux 
density. It will also affect the maps of calibrator sources that we 
took from the literature and thus indirectly affect our estimates of 
gain corrections. We do not have information about uncertainties of 
gain curves of LBA antennas. 

Table 5 displays 12 out of 633 rows of the catalogue of correlated 
flux-density estimates. The full table is available in the electronic 
version of this paper (see Supporting Information). Columns (1) 
and (2) contain the J2000 and B1950 IAU names of a source. Col- 
umn (3) contains the observation date of the experiment. Column 
(4) contains the status of the source: Cal if it was used as an ampli- 
tude calibrator. Column (5) contains the number of RCP observa- 
tions used in processing. Columns (6), (7) and (8) contain median 
values of correlated flux densities determined in that experiment 
at baseline projection lengths 0-6 MX, 6-25 MX and 25-50 MX 
respectively. Columns (9), (10) and (11) contain the arithmetic 
mean of the correlated flux-density uncertainties that accounts for 
both thermal noise and uncertainties in the calibration. If no data 
were collected to that range of baseline projection lengths, then 
— 1.0 is used as a substitute. Column (12) contains the experiment 
name. 
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Correlated flux density (Jy) 


Figure 10. The probability density histogram of correlated flux densities 
in LCS1 experiments. The dotted line with discs denotes correlated flux 
densities at baseline projection lengths 0-6 MX and the solid line with discs 
denotes the correlated flux density at projection lengths 25-50 MX, 

Fig. 10 shows the probability density histogram of correlated flux 
densities in LCS1 experiments. 


6 THE LCS1 CATALOGUE 

Table 6 displays 12 out of 410 rows of the LCS1 catalogue of 
source positions. The full table is available in the electronic version 
of this paper (see Supporting Information). Columns (1) and (2) 
contain the J2000 and B1950 IAU names of a source, respectively. 
Columns (3), (4) and (5) contain hours, minutes and seconds of 
the right ascension, respectively. Columns (6), (7) and (8) contain 
degrees, minutes and arcseconds of the declination, respectively. 
Columns (9) and (10) contain inflated position uncertainties in the 
right ascension (without the multiplier cos 5) and declination in 
milliarc seconds, respectively. Column (11) lists the correlation co- 
efficient between the right ascension and declination and column 
(12) contains the total number of observations used in the position 
data analysis, including RCP and LCP data. Columns (13), (14) and 
(15) contain the median estimates of the correlated flux density in 
Jansky over all experiments of the source in three ranges of baseline 
projection lengths: 0-6 MX , 6-25 MX and 25-50 MX, respectively. 


o 



Position error (mas) 


Figure 11. The histogram of the semimajor axes of inflated position error 
ellipses among 410 sources in the LCS1 catalogue. 

Columns (16), (17) and (18) contain estimates of the correlated 
flux-density uncertainties in Jansky. 

Of 421 sources observed, 410 were detected in three or more 
scans. In four LCS1 observing sessions, 17 731 observations out of 
19494 were used in the LSQ solution together with 7.56 million 
other VLBI observations. The semimajor error ellipse of inflated 
position uncertainties varies in the range 1.4-16.8 mas with the 
median value of 2.6 mas. The distribution of sources on the sky is 
presented in Fig. 9. The histogram of position errors is shown in 
Fig. 11. 


7 SUMMARY 

The absolute astrometry LBA observations turned out highly suc- 
cessful. The overall detection rate was 97 per cent - the highest rate 
ever achieved in a VLBI survey. If we exclude extended sources, 
non-AT20G sources and the six planetary nebulae included in the 
candidate list due to an oversight, the detection rate is 99.8 per 
cent! We attribute this high detection rate to two factors. First, 
the AT20G catalogue is highly reliable and is biased towards very 
compact objects. Selecting candidates based on the simultaneous 
AT20G spectral index proves to be a good methodology. Secondly, 
the LBA has very high sensitivity. The baseline detection limit over 


Table 6. The first 12 rows of the LCS1 catalogue of source positions of 410 target sources. The table columns are explained in the text. The full table is 
available in the electronic version of this paper (see Supporting Information). 


J2000 name 

B1950 name 

Right ascension 

^hmsj 

declination 

r „) 

cr(a) 

(9) 

(mas) 

<t(<5) 

(10) 

(mas) 

Correlation 

coefficient 

(11) 

# Observation 
(12) 

Correlated flux 
(13) (14) (15) 

(Jy) (Jy) (Jy) 

Correlated flux 
uncertainty 
(16) (17) (18) 

(Jy) (Jy) (Jy) 

J001 1—8443 

0009-850 

00 11 45.90267 

-84 43 20.0096 

46.1 

3.1 

-0.228 

30 

0.162 

0.165 

0.180 

0.005 

0.005 

0.006 

J0028— 7045 

0026-710 

00 28 41.56281 

-70 45 15.9267 

7.1 

1.7 

0.006 

40 

0.092 

0.086 

0.081 

0.005 

0.006 

0.005 

J0030— 4224 

0027-426 

00 30 17.49264 

-42 24 46.4827 

2.0 

1.4 

-0.036 

37 

0.312 

0.290 

0.256 

0.011 

0.010 

0.011 

J0033— 4236 

0031-428 

00 33 47.94145 

-42 36 14.0676 

3.7 

2.2 

-0.012 

37 

0.132 

0.135 

0.143 

0.005 

0.005 

0.008 

J0034— 4116 

0031-415 

00 34 04.40893 

-41 16 19.4729 

7.1 

3.6 

0.089 

25 

0.053 

0.045 

0.039 

0.003 

0.002 

0.004 

J0040— 4253 

0038-431 

00 40 32.51473 

-42 53 11.3916 

4.5 

2.8 

0.057 

35 

0.068 

0.058 

0.067 

0.003 

0.002 

0.007 

J0042— 4030 

0039-407 

00 42 01.22481 

-40 30 39.7419 

3.4 

2.0 

0.013 

33 

0.255 

0.253 

0.276 

0.010 

0.011 

0.013 

J0042— 4333 

0040-438 

00 42 24.86725 

-43 33 39.8164 

3.9 

2.3 

0.002 

37 

0.165 

0.113 

0.126 

0.005 

0.005 

0.008 

J0044— 8422 

0044-846 

00 44 26.68719 

-84 22 39.9895 

44.0 

2.9 

-0.237 

30 

0.278 

0.255 

0.234 

0.007 

0.006 

0.007 

J0047— 7530 

0046-757 

00 47 40.81228 

-75 30 11.3640 

6.4 

1.4 

-0.631 

41 

0.139 

0.118 

0.078 

0.006 

0.008 

0.006 

J0049— 4457 

0046-452 

00 49 16.62412 

-44 57 11.1658 

3.7 

2.1 

-0.005 

37 

0.244 

0.215 

0.207 

0.009 

0.006 

0.011 

J0054— 7534 

0052-758 

00 54 05.81337 

-75 34 03.6325 

6.3 

1.4 

-0.803 

44 

0.094 

0.087 

0.080 

0.004 

0.007 

0.006 


© 2011 The Authors, MNRAS 414, 2528-2539 
Monthly Notices of the Royal Astronomical Society © 2011 RAS 


The LBA Calibrator Survey - LCS1 2539 


2 min of the integration time varied from 7 to 30 mJy, with 7-12 mJy 
at baselines with the PARKES. 

We have successfully resolved group-delay ambiguities with 
spacing 3.9 ns for all observations. This became possible using 
the innovative algorithm exploiting the relatively low level of in- 
strumental group delay errors. 

We have determined positions of 410 target sources never be- 
fore observed using YLBI, with a median uncertainty of 2.6 mas. 
Error analysis showed a moderate contribution of the mismodelled 
ionosphere path delay to the overall error budget. Both random 
and systematic errors are accounted for in the uncertainties as- 
cribed to source positions, exploiting the overlap of 1 1 1 additional 
sources observed in LCS1 experiments with their positions known 
from prior observations. The positional accuracy of the LCS1 cat- 
alogue is a factor of 350 greater than the positional accuracy of 
the AT20G catalogue, corresponding to the ratio of the maximum 
baseline lengths of the LBA and the ATCA. The new catalogue has 
increased the number of sources with declinations < —40° from 
158 to 568, that is, by a factor of 3.5. 

We determined correlated flux densities for 410 target and 111 
calibrator sources, and presented their median values in three ranges 
of baseline projection lengths. The correlated flux density of the 
target sources varied from 0.02 to 2.5 Jy and was in the range 80- 
300 mJy for 70 per cent of the sources. The uncertainties of the 
correlated flux densities are estimated to be typically 5-8 per cent. 

This observing programme is continuing. By 2010 November, 
four additional 24-h experiments had been observed with several 
more observing sessions planned. 
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SUPPORTING INFORMATION 

Additional Supporting Information may be found in the online ver- 
sion of this article: 

Table 3. The list of target sources that were not detected in LCS 1 
VLBI observations. 

Table 5. The catalogue of source correlated flux-density estimates 
for both 410 target sources and 111 calibrator sources in each ex- 
periment. 

Table 6. The LCS1 catalogue of source positions of 410 target 
sources. 

Please note: Wiley-Blackwell are not responsible for the content or 
functionality of any supporting materials supplied by the authors. 
Any queries (other than missing material) should be directed to the 
corresponding author for the article. 
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